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Starting from the Vlasov-Maxwell system, an exact relativistic hydrodynamic closure for a special type water 
bag distributions satisfying the Vlasov equation has been derived. It has been shown that the hydrodynamic 
equations are fully equivalent to the original Vlasov-Maxwell equations. The linear stability of the exact 
hydrodynamic closure has been studied as well. It has been shown that all basic features of the small signal 
gain can be derived from the fluid dynamic description developed here. Satisfactory agreement with previously 
reported results has been also found. 
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I. INTRODUCTION 

Conventional lasers are ubiquitous sources of coherent 
electromagnetic radiation over a wide range of the spec- 
trum, from the infrared (around 10 mm) down to the ul- 
traviolet (around 200 nm). However, in the far-infrared, 
the so-called terahertz part of the spectrum (in the range 
from 30 mm to 1 mm) , or at X-ray wavelengths (less than 
10 nm), traditional lasers are not easily achievable. There 
is therefore ever growing interest in alternate sources of 
intense, ultra fast, coherent radiation covering these two 
portions of the electromagnetic spectrum. 

The free electron laseri^, (FEL) is a relatively new kind 
of laser, in which the electrons are not bound in atomic or 
molecular systems, neither they are confined by a lattice. 
The source of these free electrons is an electron beam ac- 
celerated to suitable energy in a linac or a synchrotron. 
Generally speaking, the free electron laser converts part 
of the kinetic energy of the electrons into coherent elec- 
tromagnetic radiation. This conversion is facilitated by 
a static magnetic field produced by a magnetic device 
called an undulator"^. Free electron lasers can produce 
continuous and widely tunable coherent radiation in any 
part of the electromagnetic spectrum. In addition, very 
short pulses can be emitted at any wavelength with the 
laser intensity being very high. 

A number of theoretical models have been developed 
over the years in order to describe the free electron laser. 
Although the first analyses involved techniques from 
quantum mechanics'* and quantum electrodynamics^'^, 
classical methods have proved to provide a clear and ac- 
curate picture of the underlying processes^"— . The self 
consistent Maxwell-Boltzmann equations have been fur- 
ther transformed into quasi-Bloch equations in order to 
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stress on the laser physics perspective of the probleniii. 
The plasma physics nature of the electron dynamics mov- 
ing in the undulator field has been pursued in the kinetic 
description based on the coupled Vlasov-Maxwell equa- 
tions; an approach that has been developed by Davidson 
and coworkersi^ii^. Many features of the dynamics of free 
electron lasers can be understood using a simple pendu- 
lum model. It comes about because the electromagnetic 
wave of the radiated field and the magnetic field of the 
undulator act in tandem on the electron to produce a 
sinusoidal potential similar to that of a pendulunJ^^— . 
This model is to a large extent similar to the so-called 
single-wave model, widely used in conventional plasma 
physics"'^^. 

A relativistic cold-fluid dynamics model of an elec- 
tron beam with uniform cross section propagating axially 
through a constant-amplitude helical wiggler magnetic 
field has been developed as welU^iii^. The present paper 
can be considered as a further extension in this direction. 

Starting from very basic principles of propagation of a 
relativistic electron in the static magnetic field of an un- 
dulator, and assuming that its dynamics depends on the 
longitudinal coordinate and the time only, we arrive at a 
coupled Vlasov-Maxwell system of equations in one spa- 
tial dimension and time. Taking into account an exact 
solution of the Vlasov equation in the form of a relativis- 
tic water bag distribution, we preform in Section IIIII a 
reduction of the Vlasov-Maxwell system to an exact clo- 
sure of relativistic fluid dynamic equations coupled with 
the wave equations for the radiation field. An interesting 
feature of the hydrodynamic system thus derived is the 
fact that the ponderomotive potential together with the 
pressure law enter the picture in the form of an effective 
enthalpy. In Section llVi we present the linear stability 
analysis in the derived fluid dynamic framework, which is 
fully equivalent to the kinetic description in terms of the 
self-consistent Vlasov-Maxwell equations. The numerical 
solution of the dispersion equation derived in Section IIVI 
is presented in Section IVl and certain useful properties of 
the imaginary part of its roots (the so-called small-signal 
gain) are discussed. Finally, in Section IVIl we draw some 
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conclusions. 



II. GENERALITIES AND FIRST PRINCIPLES 

The present analysis assumes a relativistic electron 
beam propagating in longitudinal s direction through a 
helical wiggler magnetic field described by 



B„ 



Bo{ex cos fco* + Gy sin kos) 



(1) 



where Bq = const is the field amplitude, Aq — 27r/fco is 
the wavelength, and and ey are unit Cartesian vectors 
in the plane perpendicular to the propagation direction. 
The vector potential associated with the wiggler field ([1]) 
is given by 



Bo 
fco 



(e^; COS kos + ey sin k^s) 



(2) 



In what follows we assume that spatial variations are 
one-dimensional in nature, so that d/dx = d/dy = and 
d/ds is generally nonzero. 

The Hamiltonian of an electron moving in the wiggler 
field, as well as in the self-consistent electromagnetic field 
can be written as 



(P + eA)^ - e$. 



(3) 



where m is the electron rest mass, 7 is the relativistic 
Lorentz factor, c is the velocity of light in vacuo and P 
is the canonical momentum. In addition, $ is the scalar 
self-potential and A is the total vector potential A = 
Aw + Ar, where A^ describes the radiation self- field. It 
is convenient to pass to dimensionless momenta and new 
" time" , dimensionless Hamiltonian, and electromagnetic 
potentials 



P 

mc 



T — Ct, 



J? = 7= Vl + (P + a) 



eA 



(4) 



(5) 



(6) 



In order to eliminate the longitudinal component of the 
vector potential as under the square root, we redefine 
the longitudinal component of the particle momentum 
by means of a canonical transformation specified by the 
generating function 

F2 (x, p; r) ^ xpj; + ypy + sps - j dsa^ (s; r) . (7) 

Dropping the tilde in what follows, we write the new 
Hamiltonian 

7 = \/l+rf + (P + a)i-(^-^y dsa,(s;T). (8) 



Here the subscript "_L" corresponds to the transverse 
components of the canonical coordinates and fields. 
Hamilton's equations of motion can be written as 



dx^ 
"dT ' 

ds 
d7 



dp J 



7c 



dr 



0, 



Ic 



dps 
dr 



ds 



(9) 



(10) 



where 7c denotes the kinetic term (the square root) in 
the total Hamiltonian ([8]). In addition, T is the electric 
force 



T ■ 



dip 



das 
97 



(11) 



acting on the particle. In the present geometry, there are 
two exact single-particle invariants in the combined ex- 
ternal A„ and self-fields A^ and $ configuration. These 
are the canonical momenta p^, transverse to the beam 
propagation direction. 

The nonlinear Vlasov equation for the distribution 
function /(x, p;t) can be written as 



a/ 

dr 



df . P±dl^f^_djAdl_ 
7c ds \ ds J dps 



7c dx±_ 
It possesses an exact solution of the form 

/(x,p;t) = S{px)S{py)F{s,ps;T), 



= 0. 

(12) 



(13) 



where p^ and Py are the exact invariants, defined by 
the last two of Eqs. ©. Expression implies that 
the transverse motion of the electrons is " cold" since the 
transverse beam emittance has been neglected. The evo- 
lution of the yet unknown function F{s,ps', t) of the lon- 
gitudinal canonical coordinates and time is governed by 
the one-dimensional Vlasov equation 



dF 

a7 



p^dF 
7c ds 



T - 



djc\ dF 



ds J dps 



= 0, 



where 7c defined as 



7c(s,Ps;t) = ^1+pI +a^(s;T), 
is the dimensionless kinetic energy for p^ 



Py 



(14) 



(15) 
with 



III. DESCRIPTION OF ONE-DIMENSIONAL BEAM 
PROPAGATION FOR A UNIFORM PHASE-SPACE 
DENSITY 

It can be verified that there exist a class of exact solu- 
tions to the one-dimensional Vlasov equation (|14|) of the 
form~i^ 



F{s,Ps;t)=c[q 



Ps ~pi Hs;t) 



e 



Ps - p. 



(+) 



(s;t)]}, 

(16) 
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where 0{z) is the well-known Heaviside step function. 
Expression (|16p implies that the distribution function has 
a constant phase-space density defined by the constant C 
within a simply connected region confined by the curves 
pi ' (s; r) and pi"'' •* (s; r) and is zero outside. If J^(s,Ps;t) 



satisfies Eq. (|T6|) initially at r = 0, then the nonlin- 
ear Vlasov equation ((T?l) assures that the phase-space 
density remains constant^^ at subsequent values of r as 
the boundary curves pi •* (s; r) and pi^-* (s; r) distort and 
evolve nonlinearly in response to the applied external and 
the self-generated fields. Clearly, the area number den- 
sity 



J dsdpsF{s,ps;T) 



(17) 



is preserved. 

Next, we derive the evolution equations for the bound- 
ary curves pi \s; r) and pi'^\s; r). Substituting the ex- 
plicit form of the solution (IT51) into the Vlasov equation 
(ITil) . we obtain 



dp'. 



(-) 



El 

7c 



^S[p 

ds 
= J" 



dT 

dpV 
ds 



6 (ps ~ p. 



(+) 



dp. 



dT 



) ds 



(+) 



-S(ps^pi-^)+S{ps-pi+^) 



-5{p,^pi-^)+6{ps-pi+^)]. (18) 



Multiplying the above equation by 1, Ps and p^ and in- 
tegrating over Ps , we find the evolution equations for the 
boundary curves in the form 



drK^' 



2 dT 



-f 

ds\ 



0, 



(19) 



dj, 



(+) 



ds 



.(-) 



dj, 



(-) 



ds 



{p^^'-pi-^y, 



(20) 



ld_ 
3d^ 



P. 



(+)2^_p(-)2 57, 



(-) 



ds 



ds 



Let us now introduce the hydrodynamic variables n, V 
and r as 



(23) 



)=|(7^^'-7^)^)=nF^2. (24) 



Here n(s; r) is the number density, V^(s; r) is the macro- 
scopic flow velocity and r(s; r) is yet unknown function 
to be specified later. From Eqs. and we readily 
obtain 

ri- = r(^ + ^), rf- = r(v'-^), ,25) 
,H..r(i.f), .■-.^r(i-f), ,20, 

Note from the above equations that the flow velocity V 
and the function F are given by 



VV = 



\(p't'^pV). r=i(7W+7r)). (27) 



From the second of Eqs. (I?f|) with due account of Eqs. 
(Uni) and (1^ . we can express the function F as 



F = 



H-a2 



where 



(l-F2)(l-24n2)' 



1 
8C2 



(28) 



(29) 



is the thermal speed squared. 

To complete the macroscopic fluid description, we need 
to express the source terms entering the correspond- 
ing wave equations for the electromagnetic potentials as 
functions of n, V and F. In the Lorentz gauge 



da^ dip 
ds dT 



= 0, 



(30) 



the scalar potential if and the vector potential a satisfy 
the wave equations 



U^^ J 



dpsF{s,Ps;T), 



(31) 



□a I 



lioe^ f dps , . 

a± / F{s,Ps;t) 

J 7c 



(21) 

where in accordance with Eq. (jl5p . we have defined 



7f ) {s;t) = ^Jl+pi^\s;Tf + a^s;T). (22) 



H ^— ^(e^ cos fcos + Gj, sin fcos), (32) 



□a. 



Poe' 



■ / dps—F{s,ps;T), 

J Ic 



(33) 
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where □ = 9^—9^ is the well-known d'Alembert operator 
and 



eBn 



(34) 



is the electron cyclotron frequency associated with the 
amplitude of the wiggier field. The integral on the right- 
hand-side of Eq. ([32]) can be expressed as 



/ 



Cln 



(35) 



while the integrals on the right-hand-side of Eqs. ((3T|) 
and p3p are standard and yield simply nT and nVT, 
respectively. 

Expressing Eqs. and (I^Hl) for the moments in 

terms of the hydrodynamic variables, we have 



|-(nr)-f |-(nn/)=o, 

OT OS 



d , , dr 
—{vr) + — ^T, 

OT OS 



(36) 



(37) 



and from Eq. ()21|) we obtain yet a third equation 



OS 



4C2 ds 



(nTV) = 2nT^VT, (38) 



which is a direct consequence of the first two Eqs. ((36|) 
and (|37p . The latter implies that the hierarchy of macro- 
scopic fluid equations is closed and Eqs. (|36|) and (|37)) 
comprise a complete hydrodynamic closure, fully equiv- 
alent to the one-dimensional Vlasov equation ([H)) . 

The macroscopic fluid equations (l36l) and (|37)) must 
be supplemented with the equations for the self-fields. 
The wave equations for the scalar potential f and for the 
longitudinal component of the vector potential Ug can be 
written in a straightforward manner to give 



Dip = 



tTiT, 



eomc^ 



Lias = nVl . 

m 



(39) 



(40) 



Note that the integral in Eq. ([35]) can be expressed as 
J ^F{s,Ps;r)^n(^l+^-vW) +0(^^y (41) 

Standard procedure is to formally apply the thermody- 
namic limit in which C grows infinitely together with the 



volume occupied by the electron beam such that the num- 
ber density n remains finite. In this approximation the 
thermal speed vt becomes infinitely small as well. To 
second order, we can rewrite Eq. ([5^ as 

m V 3 / 



ujcko 



(e^ cos kgs + ey sin kQs). 



(42) 



In order to eliminate the explicit dependence of the right- 
hand-side of Eq. on the longitudinal coordinate s, 
we introduce the helical field variables according to 

•A.X = cos kos + ay sin feps, 



Ay = ~ax sin kciS + ay cos fcos. 



(43) 
(44) 

Note that ^ al + = Al + A^ Finally, we collect 
the macroscopic fluid equations P6|) and (157)) together 
with the equations for the electromagnetic fields to write 
the basic system 



|:(,>r) + ^(„n/) = o, 



^loe'^nAx 



(45) 
(46) 
(47) 



dA 

□A-2fco V^-fcoA = 

as m 



2 2\ , ^cko 



1 H — vi^n 

3 / c 



(48) 



aAy + 2ko- 



dAx 



fcoA = ^^^^^(l + ^4«')- (49) 



Ua, = nVl . 



(50) 



of equations, which will be the starting point for the sub- 
sequent analysis. 

Worthwhile to mention is that, instead of the Lorentz 
gauge adopted in the present paper, it is possible to use 
the Coulomb gauge, where dgag = 0. In this case it is 
necessary to take Eq. (|47p in the form 



while Eq. ([50]) transforms as 



r"-r. 



nvr. 



(51) 



(52) 



It is important to note that the exact hydrodynamic 
model derived here is invariant under Lorentz transfor- 
mation. Detailed proof of this assertion is given in Ap- 
pendix |X1 
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IV. LINEAR STABILITY OF THE RELATIVISTIC 
HYDRODYNAMIC MODEL 



^ods{dr + vods)vi + -^dsivodr + ds)Axl 



Let us consider in what follows the simplest case of 
a cold electron beam, that is the limiting case, where 
nvT — >■ 0. First of all, we note that the hydrodynamic 
equations (j45|) and (|46l) admit a stationary solution of 
the form n ~ Uf) — const and V = vo — const. In addi- 
tion, the stationary transverse components of the vector 
potential are given by the expressions 



■AxQ 



Ayo = 0, (53) 



c2 1^ i?2 -^-1 



DAxi - 2kQdsAyi - klAxi = 



ni 
no 



iO^AxO 



-ni, 



DAyi + 2kodsAxi - kfAy = 0, 



where 



(61) 

(62) 
(63) 



where cup is the plasma frequency expressed as 



i.2 - t.2 , ^ 

Hi — -t- . 



(64) 



2 _ e no 



(54) 



Note also that in obtaining the stationary solutions 
above, the explicit assumption that the beam density 
no is small has been made. Therefore, the longitudinal 
space-charge effects can be neglected {(po , Uso and Tq are 
all zero). Finally, 



To = 70^0, 7o 



1 



Ro = \l + Aio- 



(55) 

Instead of the full system of equations (j^H) - (|50l) , we 
shall use 



|3(„r) + |(„rF,.o. 



drds 



[VT)- 



eo'mc^ 



(56) 



(57) 



coupled with the equations for the transverse vector po- 
tential p8)) and p9)) . Equation (l57l) has been obtained 
by differentiating Eq. (|46p with respect to the longitu- 
dinal coordinate s and making use of the Lorentz gauge 
condition (pO| and the equation for the scalar potential 

(HZD. 

Following the standard procedure, we take 



It is convenient to perform a Lorentz transformation by 
introducing a new time and longitudinal coordinate vari- 
ables according to the expressions 

^ "fo[T - vo-s), z = 7o(s-wot)- (65) 

Derivatives transform according to 

9t = 7o(90 - wo^z), ds ^ ^io{-~vode + d:,). (66) 



Differential operators entering Eqs. ([60]) and ((6T|) sim- 
plify considerably 



dr + vods = -de, 
70 



Vodr + ds 



1 

70 



(67) 



while the d'Alembert operator □ is invariant in the iner- 
tial frame moving in the longitudinal direction at speed 
Vo- Thus, in the new Lorentz coordinate system the lin- 
earized equations simplify considerably 



dgni + noJoOzVi H -^^OeAxi 



Rl 



0, 



(68) 



"fode{-vode + dz)vi 



Axo 
Rl 



dzi-vode + dz)Axi 



n2"-^2:l + V070^^1 + 
-tin 



ni 

no 



(69) 



n — no + eni + 



Ax = Axo + ^Axi + ■■■ 



V — Vo + tvi + . 



Ay f-Ayl 



(58) 



(59) 



and expand the quantities F, TV, nV and nW to first 
order in the formal expansion parameter e. Substituting 
all of the above into Eqs. ([56]), (HZ]), gH]) and dH]), and 
retaining linear terms, we obtain 

(dr + vods)ni + no^livodr + ds)vi 



»0-4a:0 ^ Vods)Axl = 0, 

Ho 



(60) 



(□ - kl)Axi - 2koio{-vod0 + dz)Ayi 



i^pAxo 
c^uo 



-ni, 



(70) 



(□ - kl)Ayi + 2fco7o(-«o5e + 9,)Ai = 0. (71) 

Manipulating the first two Eqs. (|68p and (|69|) . we can 
eliminate the linear velocity vi and obtain a single equa- 
tion relating the linear density ni and the horizontal com- 
ponent of the vector potential. We have 



{-Vods + dz) 



= 0. 



(72) 



Hydrodynamic Approach to the Free Electron Laser Instability 



6 



It is now a simple matter to obtain a single equation for 
the horizontal component of the vector potential. The 
result is 



where the operator H is given by 



(73) 



0.1 
0.08 
0.06 
- 0.04 
0.02 



Y=2, mJ(kgC)=Q.467, (mj(k^c)) =0.404, 2n/kg=l cm 




(74) 



Note that, the left- hand- side of Eq. (fTB]) represents a 
product of two operators. In the laboratory frame the 
action of the first operator {—vode + dz) simply implies 
that the linear solution Axi to our fluid dynamic model 
does not depend on the longitudinal variable s and there- 
fore is an arbitrary function of time to this end. This ob- 
servation is consistent with Eqs. (l60l) and (|6T|) . provided 
the condition 



■^xo . , 2 , '^1 n 

ST-^xi + vq^qVi H =0, 

Hq no 



(75) 



holds for the time dependent first order quantities. Note 
also that, if the velocity vi(t) is chosen arbitrarily (but 
dependent on time r only), then it-i^t) and Axi{t) can be 
determined uniquely from Eqs. (|62|) and ()75p . Namely, 
expressing ni in terms of Axi and vi from Eq. (j75p and 
substituting the result into Eq. ((62|) . we obtain 



d^Axi +ujIAxi = -^AxoVo^qVi, 



where 



2 ; 2 



(76) 



(77) 



is the characteristic frequency of harmonic oscillations. 

In the next Section, we analyze in detail the Fourier 
spectrum of the operator H. 



V. NUMERICAL RESULTS FOR THE SMALL-SIGNAL 
GAIN 



Let us rewrite the operator defined in Eq. ((74)) in the 
laboratory frame as follows 



{n-klY + Akldl 




(78) 



4 



FIG. 1. Free electron laser instability growth rate as a func- 
tion of the normalized radiation wave number. 



Next, we seek a solution to the equation "^Axi = in the 
standard form 

oo 

Axi — "^^ j dA:yl(fc) exp [ifes — iw„(fc)r], (79) 



where uJn{k) are all possible solutions of the dispersion 
equation 



-T - loikvo - cj) 



LO^-ik + kof 



uP'A^ 

^p-^xO ( 2 



up — [k — koY 




(80) 



In order to assess the predictions of the hydrodynamic 
model on a qualitative and quantitative level, Eq. (j80p 
has been solved numerically and the results are presented 
in Figures [T]- [5] Note, that in the notation adopted here 
the frequency [lo] = is, in fact, the radiation wave 
number. 

At low 7o values (Figures [T] and [2]) the observed be- 
havior of the free electron laser instability growth rate 
is qualitatively similar to that reported earlier^'^. Note, 
however, that a hydrodynamic approach is adopted here, 
whereas a kinetic description has been used by Davidson 
and Uhroi^. In addition there is a gamma factor dif- 
ference in the definitions of ujp and Wc in the paper by 
Davidson and Uhroi^ and here. 

At higher 70 and undulator magnetic field strength val- 
ues (see Figure[3]) the instability bandwidth increases and 
the location of its peak shifts towards shorter radiation 
wavelengths. 
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y=2, (0^/(/t|jf)=2.334, {aJ{k^c)T=Q.5i3, 2%/k^=l cm 



mJ(k^c)=2S, (aiJ(k^c))=l6Q, 2nlk^=2() cm 





1500 



FIG. 2. The same as in Fig. [T]but at higher magnetic field 
and beam density values. 



FIG. 5. The same as in Fig. |4]but the undulator period length 
has been increased by a factor of 2. 



7^10, mj(k^c)=l, {mJ{k^c)) =QAM, 2K/k^=l cm 



S 0.06 
E 




200 400 600 800 1000 1200 



FIG. 3. The same as Fig. [T]but at higher beam energy and 
magnetic field strength. 



(aJ(k^c)=U, ((aj(k^c)f =40, 2Tilk^= 10 cm 
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FIG. 4. Dependence of the instability growth rate on beam 
energy for fixed magnetic field, beam density and undulator 
period values. 



Figures m and [S] show the effect of beam energy on the 
instability at two different undulator period lengths. As 
can be seen at lower beam energies a single instability 
band exists, whereas at higher beam energies two sep- 
arate bands are generated. In addition (see Figure [SJ , 
there is a threshold beam energy below which no insta- 
bility exists. Increasing beam energy above its threshold 
value an increase of both the bandwidth (range of de- 
values) and the peak value of a; as a function of k can be 
observed. 



VI. CONCLUDING REMARKS 

Starting from first principle an exact relativistic hy- 
drodynamic closure of equations describing the dynam- 
ics of an electron beam propagating axially in a static 
magnetic undulator field has been obtained. It has been 
shown that the hydrodynamic equations are fully equiv- 
alent to the Vlasov-Maxwell system for a special type of 
solutions of the Vlasov equation. Moreover, as expected, 
the warm (in general) fluid dynamic equations derived in 
the present paper are invariant under Lorentz transfor- 
mation. 

Another interesting feature of our hydrodynamic pic- 
ture is the underlying pressure law. The latter together 
with the ponderomotive potential (usually proportional 
to the transverse vector potential squared) represents an 
effective enthalpy of the system [see Eq. (f28| ]. Notewor- 
thy to mention is also that in the non relativistic limit our 
system of hydrodynamic equations reduces to the well- 
known picture with triple adiabatic pressure Iaw2£i2i. 

As a direct application of the theory developed here, 
the linear stability of the exact hydrodynamic closure has 
been studied. It has been shown that all basic features 
of the small signal gain can be derived from the fluid dy- 
namic description developed in the present paper. Sat- 
isfactory agreement with previously reported results has 
been also found. 
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A possible extension of the approach initiated here 
could be a numerical modeling of the hydrodynamic 
equations, as well as analysis of nonlinear effects and 
possible formation of solitary wave patterns and coherent 
structures. All of the above we plan to perform in the 
near future. 



Appendix A: Lorentz Invariance of the Hydrodynamic 
Mode\ 



Simple rearrangement of terms in the above equation 
leads to Eq. (|l5l) in the new coordinate system. In a 
similar way it can be verified that the left-hand-side of 
the momentum balance equation (|46p does not change 
in the moving system. Since the longitudinal force T is 
an obvious Lorentz invariant due to the transformation 
law of the electromagnetic potentials (|A4p . Eq. (|46|) also 
remains unchanged in the new coordinate system. 



Let us first represent the quantity F defined by Eq. 
281 as 



F = 7S, 



where 



1 + 



1 - 7i2 ■ 



(Al) 



(A2) 



In order to prove the relativistic invariance of our basic 
system of equations psj) - (|50p . we recall the Lorentz 
velocity addition law 

^^TT^' 7 = 7o7'(l + «on, (A3) 

1 + voV 

following from the Lorentz transformation (|65p. Here, 
with prime we denote the value of the corresponding 
quantity in a coordinate system moving in the longitudi- 
nal direction s with velocity vq. In addition, the scalar 
potential f and the longitudinal component of the vector 
potential as transform according to the expressions 



a's^7o{as-vo(f), (A4) 



Using expressions (|66p for the corresponding deriva- 
tives, we obtain 

jUde-vod,)[nr'{l + voV')] 



+ j^{-vode + d,)[nr'{V' + vq)] = 0. 



(A5) 
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